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Abstract. 

We examine a sample of detected and undetected Isolated 
Neutron Stars (INSs), selected on the basis of their energy loss 
and distance, in order to maximize their detection probability. 
Since the sample we consider encompasses more upper limits 
than detections, we make use, for the first time in this field of 
astronomy, of survival analysis procedures through the ASURV 
software package (Feigelson & Nelson 1985, Isobe et al. 1986). 
We show that these techniques lead to a better understanding 
of the physical processes at work in high-energy emission from 
INSs. In particular, the recent detection of PSR 1951-1-32 and 
upper limits from pulsars pointed but not detected by EGRET 
show that the 7-ray efficiency of ISN is not correlated to any 
simple pulsar parameter. 
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1. Introduction 

After the recent PSR 1951-f32 detection (Ramanamurthy et 
al. 1995a), six INSs have been observed as high-energy 7-ray 
sources by EGRET, at energies above hundreds of MeV. One 
should also consider a preliminary report (Ramanamurthy et 
al. 1995b, these proceedings) of an EGRET detection of PSR 
0656-1-14. As discussed in a recent paper (Goldoni et al. 1995), 
before these last results, the efficiency ri-, of transforming rota- 
tional energy loss into 7-ray radiation of the first 5 high-energy 
INSs seem to grow with the decrease of the surface magnetic 
field B and of the period derivative P, and with the increase 
of the age r (Harding 1981;Thompson et al. 1994) and of the 
period P. 

As we will show, PSR 195H-32 detection rules out linear 
correlation between 77-y and some pulsar parameter. However, 
it is worth noting that the same conclusion could be achieved 
looking at EGRET upper hmits (Ulmer & Schroeder 1994), 
and even at COS-B data (Buccheri et al. 1978). 

The complete coverage of the sky by EGRET produced a 
lot of 7-ray upper limits (Thompson et al. 1994) for known 
radio pulsars in the Princeton catalog (Taylor et al. 1993). In 
order to correctly evaluate the observational data, it is nec- 
essary to take into account detections and non-detections in 



a self-consistent way (Feigelson et al. 1986): in fact the infor- 
mation involved in a non-detected source is different from the 
information coming from a non-pointed one. 

We discuss a sample of possibly detectable INSs on the 
basis of their expected energy flux at Earth E/d^ (Mereghetti 
et al. 1994, Thompson et al. 1994). This sample includes in 
the first place all EGRET pulsars with the exception of PSR 
1055-52, and all ROSAT pulsars. 

2. Selection of the sample 

A key problem is to find a statistically unbiased sample of INSs. 
This is not straightforward due to the low number of present 
day high-energy detections and to the difficulties of a complete 
radio pulsar survey (see e.g. Curtis Michel 1991). We excluded 
millisecond pulsars as no one of them was detected as high- 
energy radiation emitter (Fierro et al. 1995), neither they were 
detected as a whole in clusters (Barret et al. 1993; Michelson 
et al. 1994). 

We selected the first 40 INSs by their E/d^ ranking. This 
criterion selects young nearby objects. As it is shown in Table 
1, it works well for the first four EGRET detections and even 
for PSR 1951-1-32, but it fails to predict the detection of PSR 
1055-52. So one could expect not only high-rank INSs to be 
detectable, but even objects in a lower position (Hermsen et 
al. 1993). We took the beaming factor constant (I/Itt), consid- 
ering its variation as a second order effect (Helfand 1994). We 
excluded farther objects by simple energetic arguments, e.g. 
their 7-ray luminosity would be greater than their E with any 
reasonable beaming choice. The properties of the sources are 
listed in Table 1. 



3. Data Analysis 

An extensive field of statistics, called survival analysis, has 
been developed to deal with censored data, for which only in- 
complete information is present. It has been widely used over 
some decades in epidemiology and industrial reliability. In the 
past years several astronomers (Avni et al. 1980, Feigelson & 
Nelson 1985, Schmitt et al. 1985) applied these methods to 
samples containing upper limits. In this way it is possible to use 
each detection and upper limit in an efficient and well defined 
manner, minimizing error sources. In fact, spurious correlations 
can appear if upper limits are neglected from the analysis of 
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Table 1. Parameters of the first 40 INSs, ordered by E/d? ranking (Taylor et al. 1993). The "X" labels INSs detected by 
ROSAT. Waiting for a confirmation, for PSR 0656+14 we consider only the upper limit. Fhixcs and efficiencies are derived 
from E > 100 MeV EGRET data (Thompson et al. 1994). The photon fiux limits are 99.9 % confidence limits bsised on spatial 
analysis. The energy flux limits are obtained as in Thompson 1994. 



a sample (Elvis et al. 1981; Feigelson & Berg 1983). Simula- 
tions by Isobe et al. (1986) show that they disappear with a 
proper use of existing upper limits. For a complete description 
of survival analysis, see Feigelson et al. (1986), where a general 
bibliography as well as a discussion of astrophysical problems 
are presented. 

To perform this analysis wc used the ASURV Rev 1.2 soft- 
ware package (Isobe & Feigelson 1990; La Valley et al. 1992), 
kindly provided to us by E. B. Feigelson. The validity of these 
methods relies on random censoring: the distribution of upper 



limits must be independent from the distribution of true data 

values; in other words, high and low flux objects must have 
the same distribution of observational sensitivity. Censoring 
becomes more random if one uses not flux but intrinsic quanti- 
ties such as luminosity or spectral index, which are folded with 
other parameters like distance (Magri et al. 1988). 

However, often real observing changes all of this. When 
an observer does not detect an object, he usually lowers the 
threshold by further observations. In our case this means sum- 
ming more EGRET observation periods to achieve higher sen- 
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sitivities. The example of PSR 1951+32 detection illustrates 
this process very well, showing that three years of observa- 
tions were necessary to perform this task. We also note that 
the data reduction method of Ramanamurthy et al. (1995a) is 
different from the standard one, employing a new photon selec- 
tion method. The same can be said for PSR 0656-1-14 prelimi- 
nary detection which was obtained with a completely different 
method. 

This shows that there is usually no precise flux limit be- 
tween detections and non-detections. In the idealized case the 
sample is formed by detections above a certain level with non- 
detections below that level. In our case we will consider that, 
thanks to the great EGRET field of view, all the sky received an 
almost uniform coverage. We can so define in a satisfactory way 
an upper limit of detectable photon flux of ~ 10~^ ph/cm? /s. 
This limit is more or less coincident with the E/d? condition 
chosen for the INSs' sample. Here we use a few different tests 
to establish a correlation between INS's 7-ray emission and 
the rotational parameters. First of all we compare the result 
of simple correlation tests on the detected sources; in a second 
time we include the upper limits in the test procedure. 

It should be noted that statistical tests of survival analysis 
do not consitute a well deflned procedure in astronomical data 
analysis. This is especially true for linear fitting techniques 
which are seldom used in other fields of research. For this rea- 
son we did not follow completely La Valley et al. (1992) advice 
of using all the available tests, but we excluded Sclimitt's lin- 
ear regression method (Schmitt et al. 1985) due to unclear and 
somewhat arbitrary bin selection. We perform as a first step 
the linear correlation tests, and then apply the linear regres- 
sion tests. The second step will be mainly useful to show the 
effects of upper limit inclusion in our analysis. 



4. Results 

We first made correlation analysis on the detected sources, and 
then to the whole sample with the upper limits. We present the 
results of the same tests even on ROSAT-detccted INSs (Table 
2), and we found no difference with the results of the whole 
sample. While in the case of detected sources a correlation 
probability higher than 95% is shown by P, P and r (not by 
B) , when we consider the whole sample (or the X-ray sample) , 
only the age r remains above the 95% threshold. 

In Fig.l there are plotted the 7 efficiencies as a function of 
B and r, for the six sources detected to date and the EGRET 
upper limits for other sources in the sample. The linear fits 
for the first five detected sources (solid line) and for the whole 
sample (dashed line) are also drawn. It is apparent that the 
detection of PSR 1951+32 and the introduction of the upper 
limits greatly affect the B plot and also the r plot, reducing 
the likelihood of the fits. 

We thus conclude that a single linear fit of INS 7-ray ef- 
ficiencies as a function of any rotational parameter, with the 
possible exception of age, is not acceptable in the EGRET 
energy band. As it was shown by Goldoni et al. (1995), the 
bulk of the emission of middle-aged INSs falls in this energy 
range. The lack of detections of lower-B INSs in this band is 
a very important observational fact, as it suggests that either 
the overall electromagnetic efficiency decreases, or a great part 
of the energy loss fllls up another region of the spectrum. 
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Table 2. Correlation test results for Cox regression (1) and 
Kendall's tau (2) tests (see Feigelson et al. 1986). The quoted 
numbers represent the probability that a correlation between 
the 7-ray efficiency and the various rotational parameters does 
exist. 



Is really the overall electromagnetic efficiency decreasing? 
In the firame of the outer gap model (Cheng et al. 1986a,b) the 
answer could be yes. 

In the polar cap models (Ruderman & Sutherland 1975; 
Daugherty & Harding 1982; Arons 1983; Harding et al. 1993) 
the answer is not unique. A simple physical explanation, con- 
cerning the number of emitted photons and the cutoff energy 
of the 7-ray emission could be given. 

Old pulsars (like PSR 0950+08 and PSR 1929+10) Have 
hard spectra at high energies, with spectral indexes 5 approach- 
ing unity and possibly becoming even less than one. At the 
same time the high-energy cutoff of the emission grows, owing 
to the reduced optical depth of the magnetosphere (this is true 
if the spectral break is due to pair production rather than to 
curvature radiation cutoff). In this way we could find 7-ray 
efficiencies very high without being observed by EGRET. 

For instance, for 5 = —0.9 and Ebreak ~ 45 GeV, PSR 
1929+10 would have ij^ > 0.3 while being undetectable by 
EGRET. In fact its spectrum, combined with the efficiency, 
would give at E > 100 MeV a flux F ~ 6.6 x 10"'' ph/cn?/s, 
much lower than the Thompson et al. (1994) upper limit of 
1.63 X 10~^ ph/cm? /s. If we instead look at the flux at -B > 
1 GeV the situation is better, thanks to the lower value of 
the 7-ray background and the better instrument efficiency: the 
expected flux is F ^ 5x 10~^ ph/crn^ / s, only about four times 
less than the upper limit of 1.7 x 10~* ph/crn^ /s. The same 
arguments hold for PSR 0950+08. 

However this is not the case for PSR 1951+32, which is ex- 
pected to be a harder 7-ray emitter in this frame, with 5 ~ — 1.4 
instead of the observed 5 ~ —1.74 ± 0.11. This is a question to 
be studied separately, owing to the peculiar rotational param- 
eters of this object which led White & Stella (1988) to suggest 
it is a "recycled" pulsar. 

Clearly the way to overcome present day difficulties is to 
develop new instruments with higher effective areas in this en- 
ergy range. For example, there exists a SLAC project, called 
GLAST, which should have a 8000 cm? effective area from 0.01 
to 300 GeV, and a limiting flux Fmin = 1-5 x 10"^° ph/cm^/s 
between 1 and 300 GeV (Michelson 1995). 
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Fig. 1. 7 efficiency as a function of B and t, for tlie E/d? selected 
sample; superimposed are linear fits of the first 5 detected points 
(dashed) and of the whole sample with the upper limits (solid). 



5. Conclusions 

We investigated the high-energy 7-ray emission of the 40 INSs 
with large E/cf, all detected as radio pulsars with the excep- 
tion of Geminga. Due to the low number of detections with 
respect to non-detections, the only way to correctly deal with 
these data is to include non-detections as well as detections in 
our analysis. The presence of non-detections, in fact, rules out 
the possibility to perform a linear fit between r;^ and any INS' 
rotational parameter. 

A possible explanation of this phenomenon can be given 
in the framework of outer gap and polar cap models. In the 



first case the 7-ray emission is quenched for the majority of 
radio pulsars (Chen & Ruderman 1993), and the INS is no 
more a strong 7-ray emitter. In the second the electromagnetic 
emission, still very strong, becomes harder, emitting the bulk 
of its energy at _B > 10 GeV thus becoming unobservable by 
current instruments. Future instruments with good sensitivity 
at i? ~ 100 GeV currently being studied, will address this 
question properly. 
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